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The problem of the Kohn mode in bosonized theories of one-dimensional interacting 
fermions in the harmonic trap is investigated and a suitable modification of the 
interaction is proposed which preserves the Kohn mode. The modified theory is 



used to calculate exactly the inhomogeneous linear mobility /j,(z, Zq)u) at position z 



in response to a spatial force pulse at position zq. It is found that the inhomogeneous 
particle mobility exhibits resonances not only at the trap frequency uji but also at 
multiples me, m = 2,3, ... of a new renormalized collective mode frequency which 
depends on the strength of the interaction. In contrast, the homogeneous response 
obtained by an average over zq remains that of the non-interacting system. 
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1 Introduction 



Under suitable conditions, a gas of interacting fermions in one spatial di- 
mension can be asymptotically described by a bosonic phase operator obey- 
ing a simple equation of motion. Conditions for this correspondence between 
fermionic and bosonic descriptions are a linear dispersion of free particle states 
allowing bosonization via Kronig's identity [1,2] and the addition of an anoma- 
lous vacuum which requires large fermion numbers not to cause harm. Fur- 
thermore the bosonic Hamiltonian should be bilinear in the field operator to 
facilitate diagonalization. The best known example is the Luttinger model 
[3,4,5,6]. One-dimensional systems are interesting because even weak inter- 
actions transform the quasi-particles of Fermi liquid theory into collective 
excitations of density wave type, and correlation functions are characterized 
by non-universal power laws determined by one basic coupling constant K 
(in case of one component, i.e. spin-polarized fermions). Similar techniques 
can be applied directly to interacting Bose gases as has been done recently 
with ultracold quantum gases [7,8,9] based on pioneering work by Haldane [10] 
and reviewed in [11]. Applications to mixtures of one-dimensional bosons and 
fermions were given in [12]. Luttinger methods (for reviews cf [13,14,15,16]) 
were extended to include a trapping potential [17,18] by means of the local den- 
sity approximation. In [19,20,21,23,24] attempts were made to include exactly 
- as in the Calgero Sutherland model [25,26] - a harmonic trapping potential 
into the Luttinger approach utilizing the linear dispersion of free oscillator 
states, a concept recently adopted to quantum dissipation [27]. In the case of 
interacting quantum gases the penalty are truncated interactions in order to 
render the model exactly solvable. They then are no more faithful representa- 
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tions of real interactions. Furthermore, such type of approximation may violate 
essential symmetries of the original problem. For harmonically trapped quan- 
tum gases Kohn's theorem [28,29,30] is such a symmetry. It states that the 
response to a homogeneous external field with dipole coupling to the system 
displays a resonance exactly at the trap frequency irrespective of any trans- 
lation invariant particle-particle interactions. This raises the question which 
many-body approximations respect the Kohn theorem. This problem was in- 
tensively studied in the case of ultracold quantum gases and mainly in the 
context of bosons. In [31] it was found that the local density approximation 
at zero temperature is conserving. At finite temperature T the random phase 
approximations with exchange as well as the Hartree-Fock-Bogoliubov approx- 
imation are also conserving as argued in [32]. The problem was taken up in 
[33] and the Bogoliubov approximation was proven to be conserving at zero 
temperature. In [34] it was pointed out that the Gross-Pitaevskii equation is 
non-conserving at T > because the thermal cloud of the Bose-Einstein con- 
densate is neglected. The work [35] showed that kinetic equations extending 
the Hartree-Fock-Popov method are conserving. The general dielectric ap- 
proach was proven in [36] to be conserving too. Recently it was found [37] 
that superfluid fermions with a Feshbach resonance described by the gener- 
alized random phase approximation respect the Kohn theorem. The method 
of bosonization is intrinsically approximate and raises the same question. The 
phase formulation of a "Tomonaga-Luttinger model with harmonic confine- 
ment" was given in [24]. In studying its linear response it turned out that the 
model in its original form violates Kohn's theorem. Here, we propose a method 
to remedy this deficiency. We then go on to calculate explicitely the inhomoge- 
neous particle mobility in response to a spatially localized time varying force 
for a one-component gas. To this order we set up a linear but non-local differ- 
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ential equation for a phase operator adopted to the problem of harmonically 
trapped fermions. We solve it by means of its Green's function. 

2 The Kohn mode and bosonization 

In a strictly harmonic trap the trapping potential raises the translation mode 
at zero frequency to the trap frequency u>e independent of inter-particle inter- 
actions provided these are translation invariant. In a microwave experiment, 
when the electric field is dipole coupled to harmonically trapped charges, the 
resonance occurs exactly at this frequency. This was first noted by Kohn [28] in 
the context of magnetoresonance and later generalized to solid state systems 
[29]. It was also discussed for ultracold quantum gases [30]. In the latter case 
of neutral atomic gases the effect is described in terms of a " sloshing mode" : A 
homogeneous mechanical force alternating in time leads to a resonance exactly 
at the frequency u}£. A simple theoretical route to discuss the implications of 
the Kohn theorem is to consider the operator z s of the center of mass (cm.) 
position. In second quantization it is 



Transforming from local creation and annihilation operators to the harmonic 
oscillator representation gives 



Here, a is the inverse oscillator length and iV the fermion number. For large 
particle numbers N, a condition implicit in the method of bosonization, the 
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square roots can be replaced by yN [27]. The bosonization prescription for 
density fluctuation operators p{p) = J2 q Cg +p c q (cf. e.g. [16,15]) according to 



pip) = < 



\p\ d \p\, P<0, 
VP d+ , p > 0, 



(3) 



where the operators d and d + destroy and create collective bosonic excitations 
[6], then leads to 



1 2 } 

z s = — (di + df) = — ^ y rfMsinM0 odd (w). 



(4) 



The Fermi wave number kp is given by k F = a\ / 2N. The phase operator used 
to bosonize fermionic destruction and creation operators in the formulation 
[2] is 



u 



OO 1 

) = -i E ~r= ^ -7r<u<n, 



171=1 



(5) 



with a positive infinitesimal rj. In [23] we found that the physical phase oper- 
ator determining the particle density is, however, 



1 



<t>odd{u) = - (<f>(u) + <p + (u) - <j>(-u) - 



(6) 



= E 
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e mr] sinln-u 
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The Kohn mode requires a total Hamiltonian which gives the same Heisenberg 
equation of motion for the cm. coordinate 



1 



d t z s (t) = -^2 [[ 2 s,H 



H 



-ujz s (t) 



(7) 
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as the free Hamiltonian 



oo 

H = - huj e (d+d m + d m d+). ( 8 ) 

1 m=0 



We have introduced earlier [19,21] a Luttinger-like model of interacting one- 
dimensional fermions in the harmonic trap where the interaction is given by 

1 oo 

V c = --V c Y,m(d 2 m + d 2 +). (9) 

Z m=l 

This interaction approximates p-wave scattering in the one-component gas 
[22,23] in the simplest possible way. It is obvious that the condition (7) is 
violated because the part 

v cl = -\ v c (d\ + d\ + ) = l - hcu £ v c (d\ + d\ + ) (io) 

in the interaction results in dfz s (t) = —l 2 z s {t) with 



~2 2 

6 = UJn 



(l - Vc) + 4 (11) 



for the interacting system. Quantum mechanically, the Kohn mode is a co- 
herent state which would be turned into a squeezed state with renormalized 
frequency by the interaction equation (10). This is not admissible and one 
must subtract out that part of the interaction. We are mainly interested in 
the equation of motion for <p dd because it determines the response to ex- 
ternal perturbations in the most direct way. One then must work out the 
consequences of the subtraction prescription for the phase formulation of the 
theory. In the first step one introduces the phase operator form of V c \ using 

d, = J due lu {— dJ odd (u) + n(«)} . (12) 

— TT 
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The momentum density ft is 



i 00 



f[(u) = - — Vne~ nri sin rax (d n - dfj , (13) 
with the commutator 

[i oM (u),fl(y)\ = \&(u-v) (14) 

in the reduced interval I n = (— 7r < -u < 0). This shows that 2h(f) odd and II are 
canonically conjugate on this interval. In applications the auxiliary variable u 
is related to the physical position z inside the harmonic trap by 

(Z \ 7T 
— ) — — , z = Lp cos(«o), (15) 

where L F = y/%2N /a is half the quasi-classical extension of the Fermi sea. 
The non-linear relation between the formal variable u and physical position z 
expresses the trap topology. It is also obvious that the reduced interval I n is 
sufficient to describe the physics. 

The incriminating part in the original interaction becomes 

7T 

Vex = -v c J dudu'e iu+iu ' {^dJ odd {u)d u 4oM{u') + n(«)n(«')} (16) 

— 7T 



— huJeV c J dudv! |4sinMsinM / fl(-u)II(-u / ) cos-u cosu'd u <fi odd (u)d u i(t> odd ^ . 
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The phase form of the reduced Hamiltonian (neglecting the zero-mode which 
would cancel in subsequent calculations) is 



r 1 

H r ed = ^J du 2irKTl 2 (u) + -^—^ (d u 4> odd [ 
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~V cl , (17) 
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with 



K 



l + Vc 

\l-V c 



e = he = huji 



IK 



K 2 + l 



= Twji 



1 - VI 



(18) 



The first part of the reduced Hamiltonian is diagonal with renormalized exci- 
tation energies me, m = 1,2, 3..., while the second part is a non-local but still 
bilinear correction. The same results can be obtained by a selective renormal- 
ization of the trap frequency similar to the well known procedure in coupling 
a linear bath to a harmonic oscillator: 



(19) 



3 Inhomogeneous linear mobility 



The inhomogeneous mobility is the particle current response at position z to 
a 5- function force / 

f(z,z ;t) = F(t)5(z-z ) (20) 

at a different position z Q inside the trap. Quantitatively, the inhomogeneous 
linear mobility fx(z, Zo] oj) relates the expectation value of the local particle 
current density operator j p (z) to the force by 

oo 

G P {z)) u = J die™ (j p {z)) t 

— oo 

= J dz' n(z, z'\ uj) f^z', zo) = n(z, z ; uj) F„ 
The Hamiltonian for a time- dependent, spatially localized external force f(z, zq] t) 



(21) 
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affecting N fermions on the z-axis is 



N 



H ext (z ,t) = -F(t) " *>)■ 



(22) 



n=l 



In second quantized form it becomes 



oc 



H ext (z ,t) = -F(t) J dzi> + (z)?P(z) 



(23) 



For simplicity, a one-component gas of spin polarized fermions is assumed. 
Without enhancement of the interactions by Feshbach resonances [38] this is 
an admittedly academic case when regarding the role of interactions. The total 
particle density operator contains a part Sp s i ow which varies slowly in space. 
In an appropriate WKB expansion (for large N), it is given by [23] 



Again this has an analogous expression in the theory of Luttinger liquids and 
its bosonic equivalent [10,11]. However, in the present case, there is a non-linear 
relation between the argument uo of the phase operator and the actual spatial 
position z according to equation (15). There is also a rapidly varying Friedel 
part in the particle density which results from the confinement. Away from 
the classical boundaries z = ±L F it gives a small correction to the mobility if 
we average the latter over a small length L around z with Lp ^> L ^> l//cp. 
It will be neglected here. Utilizing the continuity equation leads to 
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(24) 
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(25) 



and shows that the slow current density is given by 
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jsiow{z, t) = -- d t (j) odd (u (z),t). (26) 

71 



We make the identification j p = j s ; ow which according to equation (21) leads 
to the generic form 

M") = (27) 

for the mobility. Coupling the external force to the slowly varying density 
amounts to the replacement 



i) + (z)i)(z) -»■ 5p dow (z) = - d z (j) odd (u (z)) (28) 

71 

in equation (23) and leads to the external phase Hamiltonian 



H ext = F(t)-4> odd (u (z )). (29) 

7T 



4 Equations of motion 

The equations of motion are most easily obtained via functional derivatives 
using 



2hd t ^ odd (u, t) = 6Hred , 2hd t IL(u, t) = 6Hred . (30) 

5U(u,t) 5<p odd (u,t) 



One obtains: 

o 

dt&odd — 27re fl — 4lUJ(V c sin u j du' sin u'fl(u') (31) 

— 7T 

and 
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dtfl = -^^- r d 2 u (t) odd 



2ixK 
1 



(32) 



7T 



u e V c sinu J du' sin u' <f> od d{u') - F(t) 7^ ^0(^0) - u). 



The equation of motion for <p odd closes at the level of the second time derivative: 



d't(j>odd = e 2 dl<f) od d - ufr'c sin u 



2f>2 



- / du' sin u'4>odd{u') 

7T 



(33) 



. x tK . . . . . 2a>£V^ . . 
F(t)— d(u {zo) —u) + F(t) — — sinuo(^o) sinu. 

ft 7T/1 



Here, essential use has been made of the relation 



r[K--)= 2u e V c . 



(34) 



The non-local operator results from the subtraction procedure and projects 
onto the subspace of the Kohn mode: This is best seen by considering ex- 
citations (p(u,t) = (4>odd{u))t i-e. averaging the equation of motion over a 
non-equilibrium state. For F = one gets 



d t <P = I d u (p-cu e V c smu 



7T 



du' sin u' (p(u', t) 



(35) 



The Kohn mode is a density (dipole) oscillation with <p K oc sin u (z) and 
associated real space density modulation Sp(z) oc d z ^p K oc zj \j\ — z 2 / L 2 F . 
It solves equation (35): d 2 (px = —<?¥k — u 2 V 2 (fK with the correct frequency 
because of e 2 = uj 2 (1-V c 2 ). The non-local projection operator has no influence 
on modes orthogonal to the Kohn mode. Their excitation energies scale with 
e, the renormalized frequency of the model. 
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5 Calculation of inhomogeneous linear mobility 



The eigenvalue problem associated with the phase equation 



2t72 . 



7T 



du' sin u' ip w {v!) 



(36) 



can be cast into a convenient form by introducing the function 



<pi(u) 



■ smw, 



(37) 



which is normalized in I n = (— n, 0) and satisfies the boundary conditions 
implied by 4> dd(u), namely (pi(— 7r) = = <pi(0). Introducing the linear sym- 
metric operator 



2 2 i i 2 

6 , ^ -a 



UJ 



(38) 



the eigenvalue problem becomes 



L^n = -A^(w)y n . 



(39) 



Because of the relation e 2 = uj} (l — Vf^j it is easy to see that ipi(u) is an 
eigenfunction with eigenvalue 



= 1 - - 2 . 



la-',- 



(40) 



Furthermore, 



<£n(u) = \l — sinnw, n = l,2,... 



is an orthonormal basis on It is thus clear that for n > 2 



(41) 
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/ 2 2 i , 2 \ i i 2 

L w </?n = \-n — + — </?„ = -A n (w) <p n , \ n {uj)=n — f (42) 

holds. The relevant Green's function is 



G u (u,u) = -}^ . (43) 

„=1 A nl^J 



This allows to solve for the inhomogeneous local mobility according to equation 
(27) with (<f) dd)iu = iPu>(uo(z),u (z )). Utilizing the identity IK = uj e (l + V c ) 
leads to 



F^ue fi(uo(z))ip 1 (uo(z )) _ F^eK ~ ip n (uo(z))ip n (u (zo)) 
<Pu '~ h uJ-uj 2 h ^ n 2 ~e 2 -u 2 • { ' 



Using sin-u (2;) = —\Jl — z 2 / ' L 2 F = —Z(z) finally gives 



liujuji Z(z)Z(z ) 2icueK sin nu (z) sin nu (z ) . . 



Analyticity of the response function for ^suj > requires uj — > + «?y with a 
positive infinitesimal 77. It is seen that the inhomogeneous mobility has res- 
onances at all excitation frequencies including uji. This applies also to the 
non-interacting case. It is only the homogeneous response which is trivial- 
ized by Kohn's theorem. In order to see this, we calculate the homogeneous 
mobility via 



fj,(z,u) = —L F J du sin u /i(z, z (u );uj), (46) 
and find 
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This is identical to the local response of non-interacting fermions in the har- 
monic trap to a homogeneous force. The Kohn theorem is clearly fulfilled. 

Returning to the inhomogeneous mobility, the sum can be evaluated analyti- 
cally. By including the n = 1 term in the summation one obtains 

%K 

(j,(z,zo;u) = - . — r (48) 
Z7m sm(7ra) 

x{cos [a(7r + Uo(z) + u (z ))] - cos [a(ir - \u (z) - «o(^o)|)]} 
2V C Z(z)Z(z ) f icuLu e ((l + V c )uj - uj 2 ) 1 

with a = u/e. Here, the first term also includes a contribution from the Kohn 
mode ip i which removes all renormalizations due to the interaction from the 
second term in going over to the homogeneous response. By this mechanism 
Kohn's theorem is restored. 
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